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SUMMARY 


The unsteady transonic equation for oscillating thin 
wings is solved by a direct finite difference method in the 
case where the steady flow affects only one coefficient in 
the equation. Both concave and delta wing planforms are 
solved and the program may be used for relatively arbitrary 
planforms. Both pitching and plunging modes are calculated 
for a reduced frequency range from .2 to 1.0 . 

The results are consistent with earlier asymptotic inves- 
tigations and are in numeric agreement within the order of 
accuracy of these solutions. The thickness effect while 
small, increases as reduced frequency decreases to . 2 . 
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1. INTRODUCTION 


The prediction of unsteady flow requires taking into 
account the corresponding steady flow. Landahl’s solutions 
of steady flow completely separated the two effects thus 
omitting the effect of body thickness on the unsteady behavior 
of the body. Total inclusion of the terms omitted by Landahl 
on the other hand gives a difficult equation of mixed type 
similar to the non-linear steady transonic equation which must 
usually be solved by time consuming numeric methods. The only 
simplification is that the flow-type regions are identified 
in advance. But mixed type differencing is still required and 
^^^fher complicated by the presence of in— phase and out-of- 
phase components of the solution. In addition, a multiplicity 
of solutions is needed for one body to account for several fre- 
quencies and for various modes of motion. Thus a method which 
is fast and yet provides some measure of the effect of the 
steady flow is badly needed. 

The present report describes a numerical solution procedure 
of a simplified unsteady transonic equation which is very fast 
and reasonably accurate and which takes into account many of 
the effects of the steady flow field. The numeric solution of 
this equation is accurate within a few percent and can be 
accomplished on an IBM 360/65 computer in less than one minute 
per case (one frequency and one mode of oscillation). it is more 
flexible than the older transonic box method^ in that it 
accommodates rather arbitrary planform shape and is easily 
capable of handling variable local Mach number effects from 
the steady flow. 

2. PROBLEM FORMULATION 

Consider a rigid pointed wing which performs harmonic 
pitching oscillations of small amplitude in a steady uniform 
transonic flow. It is required that the wing be smooth and 
sufficiently thin that the conditions for the validity of small 
perturbation theory are satisfied. 
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With the flow directed along the positive x axis and the 
wing oriented along the x axis, we may introduce non-dimensional 
variables x = bx, y = by, z = bz, and t = (U^/b)t where b is 
the body length, is the free stream speed and the bar coordi- 
nates are the physical coordinates. The steady state position 
of the wing, which for convenience is assumed to be symmetric, 
is then as shown in Figure 1. Denoting the thickness ratio by 
e, the requirement that the wing be thin becomes e«l. 

Since the shock strength is or order e we may ignore ' 

the irrotationality effect of the shock effect in the solution 
2 

up to order e ^ne. Consequently , we may assume the existence 
of a velocity potential $ such that the x, y and z components 
of the flow velocity are (1 + ‘J’x? » J^’espectively:’ 

It has been shown (e.g. by Landahl^^J, that the potential o 
must satisfy the following equation: 


' ' 2 ' P 

, (1 - M ) + - 2M ,4^^. 

- *tt ' 


( 1 ) 


Consistent with the requirement that the wing perform , - 
small oscillations of amplitude 6 about its steady state 
position, we may write the equation for points on the wing as 

, i IcI" ^ 

z = e g(x,y) sgn z + 6 Re (e : f(x,y)) (2) 


where g(x,y) is the steady state wing shape of order 1 and 
f(x,y), also of order 1, represents the change in shape due 
to oscillation, k is the reduced frequency of oscillation 
equal to OJb/U^ where oo is the physical frequency. We assume 
the oscillation to be a small disturbance to the steady state 
solution, so e»6 . 

The condition of tangential flow on the body becomes: 
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( 3 ) 


i 

] = €g^ sgn z + 6 Re ( (f + ikf)e ) + G(e,8) 

^ Jbody ^ ^ 

ikt i 

where G(e,6) *= 0 - (€g sgn z + 6 Re '(.f -e. . );> 

^ Jbody ^ ^ 

» 

Letting a denote the semi-span to chord ratio, we let 
y = y/a in the last term in (3). 

G(e,5) = $„ 1 (e a ^ g- sgn z + 8 a ^ Re( e^^^ f— ) ) (4) 

^ Jbody ^ y 

Then, if the wing is assumed to be nearly planar so 
that a»€ , (4) is negligible and the boundary condition (3) 

becomes 

& ] = Re {e g„ sgn z + 8 (f + ikf) e^^*} (5) 

=^Jbody ^ 

Moreover, we assume that 6 is sufficiently small that 

(5) may be evaluated on z = ± 0 instead of (2). 

The non-linear term in equation (1) is the fundamental 

difficulty in solving the above problem. Following the 

ideas of Lin, Reissner, and Tsien^-^^ , Landahl showed that the 

non-linear term would be negligible for k sufficiently high, 

-1 -1/3 

i.e. (k»ae Ina e ' ). This restricts the range of validity 

and leaves out completely the effect of body thickness as the 
equations for the unsteady part of the potential are then 
completely independent of the steady state solution. Teipel^^J 
extended an approach used by Oswatitsch and Keunel^^^ for steady 
sonic flow by approximating the o term with a parabolic 
constant r ^^transmitting the effect of the steady flow into the 
unsteady equation. This method has been exploited by Kimble, 
Liu, Ruo and WuC®^ to obtain asymptotic solutions for low 
aspect ratio pointed wings. 

Instead of this approach $ may be represented as the sum 
of a steady potential and an unsteady potential ^2 
assuming the unsteady flow is: a sii!ia:ll?'perturbation?'.of.o.,the;fiiS43eady 
flow gives, as in 
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( 6 ) 


<1 - ^2xx - «lxx»*2x + **lx»*2xx> 

■*■ **2yy '*' **2zz “ ^2xt “ ** **2tt '■,® 


In sonic flow 1 - M = 


lx 


). It has become common 


practice to ignore the term /^ix^2xx comparison to 

0 0^ . Th'isr reflects in part the vanishing of the receding 

Xxx "I 

wave as M->1. Kimble and Wul^ carried out numerical solutions 
of two dimensional wings with and without the term 
using Spreiter’s^-^^ steady state approximate solutions and 
found less than one-half percent difference in the solutions, 
well within the order of numeric approximation. This approxi- 
mation requires further study to place it on a firmer physical 
foundation. The resulting governing equation for 02 is 

[M^ 1) ^2^ = ^2yy '>2X2 ' 

Letting 02 = 5 Re (cp and denote the Laplacian with 

respect to y and z we have finally; 


1) ‘I’lxx T + 


I 2^.2 

k M 
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with boundary conditions 


9z (x,y, +0) = (f + ikf), y| <<s(x) 


(9) 


cp (x,y, + 0) = 0 , [y[ > s(x) 

9 (0,y,z) = 0, z>0 

limcp = 0, x>0. 

2 2 
y + z ^ CO 
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3. SOLUTION METHOD 


Ignoring the fact that the leading coefficient is complex, 
equation (8) resembles a heat equation with x playing the role 
of time and spatial variables y and z. (It is actually a 
biharmonic wave equation if the complex coefficient is eliminated’ 
by dif ferentat ion) . This suggests an adaiptation of the Crank- 
Nicolson method for solution. 

The boundary condition on would behaveaas if one were 
heating a plate at the edge over an expanding region. This 
causes severe difficulties in the use of a difference method 
mesh which is uniform in x. The mesh would not match the wing 
planform edges exactly (Figure.^ 2a) . Nonetheless, this approach 
was completely implemented and filled badly. Severe oscillations 
were present in the numeric solution on the wing surface. These 
oscillations which are normally a sign of instability in the 
numeric method did not extend into the field. In addition, a 
■Von-Neumann stability analysis showed the field equations were 
stable. A similar difficulty was experienced with the transonic 
box method^^^ and was partially repaired by an edge correction. 
Various corrections were tried but a difference scheme can only 
recognize a change of the boundary to an accuracy the same as 
mesh size. This approach was therefore abandoned. 

A transformation used by Landahl^-^^ to transform a delta 
wing problem gave the fundamental idea for a means to transform 
the problem in such a way that the mesh aligned itself with 
the planform edges thus eliminating the boundary instability 
problem. 

Let y = ys(x) (10) 

z = zs (x) 

Then the equation (8) becomes) 

ct [s^ (p^ - s's (y cp- + z op~)] + p s^ cp = cp (11) 
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where [M^ (7 + 1 ) + 2 M?ki] and p while 

w = f + iki 

X 

The boundary conditions (8) become 

cpg- (x,y, + 0) = ws |y| < 1 (12) 

cp (x,y, + 0) = 0 |y| ^ 1 

and the other boundary conditions are unchanged. See 
figure 2b. 

This problem was programmed using a relatively straight 
forward adaptation of the Crank-Nicolson method for variable 
coefficients and complex cp. The stability was excellent and 
the previous difficulty completely overcome . Unfortunately, 
an iterative method of solving the equations at each x step 
was used. The amplification factors found by the Von-Neumann 
analysis approached one #or x near the trailing edge and 

= 0 . The iteration procedure took more and more iterations 

XaX ^ , 

to converge as x approached 1 . 0 . Under re.laxa!tioh..i lead<.; toi.o * 
eventual convergence but only after consuming a large amount 
of computer time. The effects were not as severe with ^ 0 

but efficiency still suffered. 

Since reprogramming was again necessary to use a direct, 
elimination procedure, an additional transformation was also 
tried to reduce storage requirements. Many mesh points must 
be used off the wing in the lateral as well as the vertical 
direction to take account of the field effects. The transfer- 
matlon _ _ 

y = sin y cosh S’ 

z = cos Y sinh ’Z’ 

eliminates the need for mesh points displaced laterally off the 
wing and in addition smoothesothe edge singularity in the 
derivatives of cp giving even better numeric accuracy there. 

See figure 2 c. 
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The equation and boundary conditions (with the double 
bars dropped) become 

2 2 2 * 

a [fcos y + sinh z) s cp - ss (sin (2y) cp /2 + sinh? (2z) 9 /2)J 

X y 2 

2 2 2 2 

+ p (cos y + sinh z) s 9 = V 9 ( 14 ) 

9 z = s w cos y, |y| < | ( 15 ) 

9 (x, z) = 0 , z > 0 

cp (0,y,z) = 0 , z > 0, |y| < Z 

lim 9*=0, x>0, |yl<§ 


This method gave good stability and accuracy of the in- 
phase part of cp. However, the out-of-phase part of 9 is iOhe? 
orders) of magnitude smaller than the in-phase part. Without 
fine mesh the out-of-phase part was lost in noise. To correct 
this 9 was split into two parts, one a known function related 
to the slender body solution, and the others' of^jcomparable 
magnitude!;;/in both in and out-of -phase parts. 

■/r . 

Cp as ■• 0 '+!,^ w cos(y)e; ^ ( 16 ) 

This was substituted into ( 14 ) and ( 15 ) and only i/' was computed 
numerically resulting in very accurate solutions for both in 
and out-of-phase parts of cp. 

4 , PRESSURE COEFFICIENT AND GENERALIZED FORCES 

The pressure coefficient is computed using central differences 
to approximate the derivatives of 9 from 

= -2 (cp + ikqp) ( 17 ) 

P X 
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We define the generalized aerodynamic force coefficients 


by 



® ’ z=0 


(18) 


where S is the wing planform area, the displacement distribu- 
tion function for harmonic oscillations in mode j is given by 
f^(x,y) cos kt (f. is of order 1), and cp. is the unsteady 
potential due to mode i. The integrals were evaluated Xafter 
integration by parts to eliminate derivatives) by Simpson^s 
rule. 


5 . RESULTS 

5. 1 Comparison With Transonic Box Method 

The program was tested for a case solved by Rodemich 

r 2 1 

and Andrew*- J by the transonic box method. The planform 
was a delta wing of aspect ratio 1.5 oscillating in plunge 
at a reduced frequency k of .5 . Results are plotted in 
Figure 3 and show agreement to within about 3%. The finite 
difference solution uses 16 points in the x direction andv ! 

8 in the y direction. In a sense this is comparable to 128 
boxes in the transonic box method. The transonic box solution 
shown, however, required approximately 300 boxes on the wing. 

The far field was approximated by setting the correction 
potential to 0 at z^ approximately 4. This means that at 
z^ = 4 the potential is forced to agree with the slender body 

solution. A value of z =8 caused only about one percent 

00 

change in the solutions whereas z^ s 2.0 caused much more severe 
variation| 

Similar tests were made with varying mesh sizes. Doubling 
the number of mesh points in any direction again caused only 
about one percent variation in the solutions. 
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5.2 The Delta Wing With Thickness Effect- 


The generalized moments and Lg 2 were computed for a 
delta wing of aspect ratio 1.5 in pitching and plunging modes 
at reduced frequencies .2, .5, .8 and 1.0 with and without a 
thickness effect. 

The thickness effect chosen for this study was based on 
earlier work of Liu and Ruot®! which extended the earlier 
work of Teipel mentioned above. It was assumed that e = .0683, 

O' - .375 and that the wing was a biconvex with simple wedge 
chord sections. This implies that the equivalent body is a 
cone of base area .0683. According to [9] this gives a 
parabolic constant r = . 55 which serves as an approximation 
to the steady state term (7 + 1) M JZ^^xx* flow was 

assumed to be sonic (M = 1) everywhere for this report. 

There is no limitation on the program (or its efficiency) to 
constant values of (7 + 1) M > however. Better data can 

be incorporated in future studies. 

The results are shown in Figures 4 and 5. Note the expanded 
scales. The effect of thickness is seen to increase as reduced 
frequency decreases which is in agreement with past asymptotic 
analysis of thickness effect^^^’ The near agreement with 

and without thickness effect at k = 1.0 is rather coincidental 
although as frequency increases, thickness effect will vanish. 

The accuracy suffers somewhat near k = 1.0 because of the 
relatively large mesh size. 

The phase angle is affected only very little by thickness. 
Note that L 22 seems toushow very little thickness effect. This 
can be traced to the fact that 9 is not normalized by k as it 
is in the plunging case. We suggest using the downwash w = ikx 
and normalizing tp to accurately account for thickness effect 
in this case. 
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5.3 The Concave Wing With Thickness Effect 

To illustrate the capabilities of the method to handle 
general planforms, a concave wing with planform given by 

s (x) 0 (x + x ^)/2 

and 0 = .4167 corresponding to aspect ratio 2.0 was chosen. 

The results show a slower onset of thickness effect with 
decreasing frequency and a slower decay as frequency increases. 
The first tendency is again consistent with qualitative trends 
noted in the earlier asymptotic investigation at low frequency. 
The second tendency has not been seen before and requires 
further investigation. Again phase angle is not greatly 
affected and the pitching results are largely swamped due 
to lack of normalization. 

5.4 Pressure Coefficients 

Pressure Coefficients were calculated for both wings 
oscillating in the plunging mode at k = . 2 . The amplitudes 
of pressure show little change when thickness effect is 
included. Thickness effect is largely reflected in rather 
significant changes in the phase angles. This emphasizes the 
importance of calculating the but—of— phase part of ip correctly 
since this part almost totally determinfes the phase shift due 
tbo thickness effect. 

6 . CONCLUDING REMARKS AND RECOMMENDATIONS. 

An efficient and accurate method has been developed for 
solving a simplified version of the unsteady transonic equation. 
The equation may be modified to include the effect of 
wherever this is not negative without interfering with efficiency 
or accuracy. Such an approach is reasonable in sonic flow when 
the receding wave effects can be ignored. In consequence wake 
influence isj ignored and shocks must be weak. Nevertheless, 
when the sonic pocket extends over most of the wing, these 
assumptions are very nearly fulfilled. The present method 
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can easily account for variable local Mach number and rather 
arbitrary planform so long as the basic assumptions are 
fulfilled. 

Additional convergence studies in the case of curved 
planforms would be valuable since accuracy is somewhat affected 
by such changes. The present program may be modified to use 
approximately 70% less storage and about the same proportion 
less time to make such studies feasible. 

It is recommended that pitching results be compared using 
a downwash function w =* ikx rather than the more conventional 
w « 1 + ikx since the present procedure gives solutions 
largely dominated by the plunge-like case w = 1. Comparison 
of results calculated with w ikx will allow normalization 
by frequency k and will more accurately reflect differerences 
due to thickness, planform, and local Mach number. 
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Figure 3. 


Real And Imaginary Parts Of the Unsteady Potential 
<p In The Plunging Mode For An Aspect Ratio 1.5 
Delta Wing At k » 0.5f Chordwise Distribution For 
.25. 























250 

240 

23 . 0 . 

^22 

220 

to 

“ 21 0 
200 
190 
180 



To Pitch For A Concave Wing. 














REFERENCES 


[ 1] Landahl, M. T. , Unsteady Transonic Flow , Pergamon 
Press, New York, 1961. 

[ 2] Rodeittich, E. F. , and Andrew, L. V.* "Unsteady 

Aerodynamics for Advanced Configurations, Part II — 

A Transonic Box Method for Planar Lifting Surfaces." 
FDL-TDR-64-152, Part II, May 1965, Dayton, Ohio. 

Lin, C. C. , E. Reissner and H. S. Tsien, "On Two- 

Doraensional Non- Steady Motion of a Slender Body in 
a Compressible Fluid," J. Math. Phy. , Vol. XXVII, 
1948, pp. 220-231. 

[ 4] TeipeiL, I. , "Die Instationaren Luftkrafte bei der 

Machzahl 1," Z. Flugwiss. , 12, Heft 1 (1964); Also 
"The Nonsteady Aerodynamic Forces at Mach Number 
One," Translation, RSIC-808, Redstone Arsenal, 
Alabama. 

[ 5] Oswatitsch, K. , and F. Keune, "Flow Around Bodies of 

Revolution at Mach Number One," Paper Presented at 
Brooklyn Polytechnic, Conference on High-Speed 
Aeronautics, January 20-22, 1955. 

[ 6] Kimble, K, R. , D. D. Liu, S. Y. Ruo, and J. M. Wu, 
"Unsteady Transonic Flow Analysis for Low Aspect 
Ratios Pointed Wings", AIAA Journal, Vol. 12, No. 4, 
April, 1974, pp. 516-622. 

Hi 7J Kimble, K. R. and J. M. Wu, "An Approximate Solution of 
Unsteady Transonic Flow Problems," Technical Report 
AFFDL-TR-74-32 , To appear as An Air Force Contractor 
Report, Air Force Flight Dynamics Laboratory, Air 
Force Systems Command, Wright -Pat ter son A.F, B. , Ohio, 
April, 1974. 


27 



[ 8] Stahara, Stephen S, and John R. Spreiter, "Development 
Of a Nonlinear Unsteady Transonic Flow Theory", 

NASA Contractor Report, NASA CR-2258, June, 1973. 

[ 9] Ruo, S. Y- and D. D. Liu, "Calculation of Stability 
Derivatives For Slowly Oscillating Bodies Of 
Revolution At Mach 1.0", Contract NAS8-20082 , Prepared 
For National Aeronautics and Space Administration, 
Marshall Space Flight Center, Alabama 35812, February 
1971. 


28 



